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Abstract. 

We present two new dynamical tests of the biasing hypothesis. The 
first is based on the amplitude and the shape of the galaxy-galaxy cor- 
relation function, £ 9 (r), where r is the separation of the galaxy pair. 
The second test uses the mean relative peculiar velocity for galaxy pairs, 
v 12 (r). This quantity is a measure of the rate of growth of clustering 
and it is related to the two-point correlation function for the matter den- 
sity fluctuations, Under the assumption that galaxies trace the mass 
(£ 5 = the expected relative velocity can be calculated directly from the 
observed galaxy clustering. The above assumption can be tested by con- 
fronting the expected v 12 with direct measurements from velocity-distance 
surveys. Both our methods are checked against N-body experiments and 
then compared with the and v 12 estimated from the APM galaxy 

survey and the Mark III catalogue, respectively. Our results suggest that 
cosmological density parameter is low, Q m ~ 0.3, and that the APM 
galaxies trace the mass at separations r > 5/i _1 Mpc, where h is the 
Hubble constant in units of 100 km s _1 Mpc. The present results agree 
with earlier studies, based on comparing higher order correlations in the 
APM with weakly non-linear perturbation theory. Both approaches con- 
strain the linear bias factor to be within 20% of unity. If the existence 
of the feature we identified in the APM £ s (r) - the inflection point near 
£o = 1 - is confirmed by more accurate surveys, we may have discovered 
gravity's smoking gun: the long awaited "shoulder" in £, generated by 
gravitational dynamics and predicted by Gott and Rees 25 years ago. 



1. Introduction 

We have been recently working on two projects, directly related to the topic 
of this Meeting - the transition to nonlinearity in gravitational clustering. We 
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have therefore decided to present some of our results here. The first idea - to 
use the inflection point in the galaxy-galaxy correlation function to constrain 
biasing is in a more mature state than the second, which uses measurements of 
the relative motions in pairs of galaxies. Accordingly, the first project is almost 
a paper - it has been submitted for publication in Monthly Notices. The results 
of our work on the second project are still far from completion but ripe enough 
to be presented and discussed. 

1.1. The CDM crisis and biasing 

The concept that galaxies may not be fair tracers of the mass distribution was in- 
troduced in the early eighties, in part in response to the observation that galaxies 
of different morphological types have different spatial distributions, hence they 
cannot all trace the mass (there are two excellent reviews on the subject: Strauss 
& Willick 1995 and Hamilton 1998). However, there was also another reason: 
to "satisfy the theoretical desire for a flat universe" (Davis et al. 1985, p. 391). 
More precisely, biasing was introduced to reconcile the observations with the 
predictions of the Einstein-de Sitter cold dark matter (CDM) dominated model. 
At the time, it seemed that just a simple rescaling of the overall clustering ampli- 
tude by setting £ ff (r) = b 2 £(r), where b « 2 might do the job (Davis et al. 1985). 
However, very soon thereafter, it became clear that this is not enough: while the 
unbiased (b = 1) £(r) had too large an amplitude at small r, the biased model 
did not have enough large-scale power to explain the observed bulk motions (Vit- 
torio et al. 1987). A similar conclusion could be drawn form comparison of the 
relative amplitude of clustering on large and small scales (eg Maddox et al 1990). 
The problem with the shape of £(r) became explicit when measurements of £ 9 (r) 
showed that the optically selected galaxies follow an almost perfect power law 
over nearly three orders of magnitude in separation. This result disagrees with 
N-body simulations. The standard (f2 m = 1) CDM model as well as its various 
modifications, including VL m < 1 and a possible non-zero cosmological constant, 
fail to match the observed power law (see Fig 11-12 in Gaztahaga 1995, Jenk- 
ins et al. 1998; most of these problems were already diagnosed by Davis et al. 
1985). Two alternative ways out of this impasse were recently discussed by Rees 
(1999) and Peebles (1999). We believe that it will be helpful to discuss both 
approaches because their existence provides the motivation for our work. 

1.2. Environmental cosmology 

A possible response to the CDM crisis is to build a model where simple phe- 
nomena, like the power-law behavior of £ 9 are much more complicated than they 
seem. In particular, one can explore the possibility that the emergence of large 
scale structure is not driven by gravity alone but by "environmental cosmology" 
- a complex mixture of gravity, star formation and dissipative hydrodynamics 
(Rees 1999). A phenomenological formalism, appropriate for this approach was 
recently proposed Dekel & Lahav (1999). According to the old, "linear biasing" 
prescription, at each smoothing scale, the galaxy- and the dark matter-density 
fields, S g and 5, are related by the linear function 

S 9 = bS, (1) 
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where b is a time- and scale-independent constant. In the new picture, the rela- 
tionship between the two fields is neither linear nor deterministic. The biasing 
factor, here defined as the ratio of the rms values of the two fields, b = <7 g /a is 
allowed to depend on the smoothing scale as well as on the cosmological time. 
A convenient measure of the stochasticity in the relationship between the two 
fields is the cross-correlation coefficient, 

R = ^ ; |*| < 1 . (2) 

The special case R = 1 describes the deterministic bias, while R = corresponds 
to the "maximum stochasticity" , or the case when the two fields are completely 
uncorrelated. The parameters R and b describe he biasing and stochasticity in 
the galaxy distribution relative to the spatial mass distribution. When referring 
to specific galaxy surveys, we will sometimes use subscripts, e.g. bjRAS f° r 
the IRAS survey. These quantities should be distinguished from the bias and 
stochasticity measures for two classes of galaxies of different morphological types, 
e.g. for early (subscript e) and late (subscript I) galaxies: b e i = a e /cri, and 

Rd = (<5 e ^)(cj e cj;) _1 - 

The above parameters can be estimated or constrained by more or less di- 
rect measurements from galaxy redshift surveys, peculiar velocity data and other 
observations. They can be also studied in semi-analytic theoretical models (Sel- 
jak 2000, Scoccimarro et al. 2000), hydrodynamic simulations or semi-analytic 
models combined with N-body simulations (eg. Blanton et al. 2000, Somerville 
et al. 2001 and references therein). 

1.3. Constraints on biasing 

If biasing is indeed important and complicated, we should expect that b,b e i,R 
and R e \ are all significantly different from unity and scale-dependent. As we 
show below, at sufficiently large scales, r > 10/i _1 Mpc, there is actually evidence 
to the contrary: the admissible deviations of b and R are small and always 
comparable to the accuracy of the measurements. 

Skewness and higher moments. The strongest constraints on large 
(weakly non- linear) scales come from the measurements of the two-, three- 
and four-point connected moments of the density field in the APM catalogue. 
These constraints are obtained as follows. One assumes that galaxies trace 
the mass and the large-scale structure we observe today grew out of small- 
amplitude, Gaussian density fluctuations in an expanding, self-gravitating non- 
relativistic gas. If our assumption is correct, by now nonlinear gravitational 
instability would have driven the distribution away from gaussianity, generating 
skewness and higher connected moments, which can be calculated analytically 
(Juszkiewicz et al. 1993, Bernardeau 1994a, 1994b). The assumptions about 
initial gaussianity and lack of biasing can then be tested by comparing the an- 
alytical predictions with observations. The predictions are in good agreement 
with the data from the APM (Gaztahaga 1994, Gaztahaga & Frieman 1994, 
Frieman & Gaztahaga 1999) as well as the IRAS PSCz catalogue (Feldman et 
al. 2000) and suggest that b{r) is within 20% of unity for r > 10/i~ 1 Mpc, or at 
linear scales (where the clustering amplitude is less than unity). 
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Redshift distortions. Some of the most radical claims that b e \ can be as 
large as 1.5 to 2 are based on comparisons of the estimates of the strength of 
clustering in the 1.2 Jy IRAS catalogue with optical redshift surveys (Strauss & 
Willick 1995 and references therein). Indeed, those earlier studies, summarized 
by Hamilton (1998), gave a redshift distortion parameter Piras = ^m 6 /^iras 
0.77 ■ 0.22. while the average and the standard deviation of the same parameter, 
obtained from optical surveys is, according to the same compilation, floptical = 
0.52 ± 0.26, implying a relative bias b opt icai/biRAS = b e i « 1.5. The most 
recent analysis (Hamilton et al. 2000) of the larger PSCz sample [Saunders et al. 



2000| gives Piras = 0.41 ± 0.13, consistent with b e i = 1. Moreover, Hamilton 



et al. also conclude that their results are consistent with RiRAs( r ) ~ 1 at 
r > 10/i~ 1 Mpc. More recently Hamilton and Tegmark (2000) have studied 
the large as well as small-scale clustering in real space, reconstructed from the 
redshift space data of the PSCz survey and found evidence for scale-dependent 
bias; however, the effect is significant at small, strongly-nonlinear scales only. 

Large scale flows. Recent measurements of the mean relative pairwise 
velocity of galaxies allow an independent estimate of f2 m and the biasing param- 
eter. These results are consistent with £l m ~ 0.3 and R ~ b ~ 1, b e i ~ R e i ~ 1 
at separations r > 5/i _1 Mpc (Juszkiewicz et al. 2000a). 



Weak lensing. The values of b and S7 m , derived from large scale flows 
are consistent with recent measurements of cosmic shear correlations based on 
the VIRMOS deep imaging survey (Van Waerbeke et al. 2001). On much 
smaller scales, r ~ l/i _1 Mpc, weak lensing data from the Sloan Survey give 
Q, m R/b ~ 0.3 (Fischer et al. 1999), again in agreement with the estimates from 
the pairwise motions, although unlike the relative velocity approach, the Sloan 
results suffer from degeneracies between O m and bias. 

Qualitative arguments. Qualitatively, strong biasing effects would be 
difficult to reconcile with the well known fact that the galaxies, dwarf galax- 
ies and IRAS galaxies have strikingly similar distributions, all avoiding the voids 
(Peebles 1993, pp. 638-642). Another empirical argument against biasing is pro- 
vided by the universal character of the observed Tully-Fisher and fundamental 
plane relations (see, for example Binney 1999; Peebles 1999). 

Simulations. All simulations of the galaxy formation process, either of 
hydrodynamic or semi-analytic variety predict morphological segregation as well 
as b and R parameters dependent on scale and time. Since these models are at 
a relatively early stage of their development, the details differ from model to 
model on small scales and at high redshifts, i.e. exactly where their b and R 
parameters are significantly different from unity and where clustering is strongly 
nonlinear. However, most of these numerical experiments broadly agree that 
with decreasing redshift and increasing scales, b and R approach unity (see eg. 
Figure 18 in Somerville et al. 2001 and references therein). 

In the end, it may turn out that to explain the observed structure, all 
we need is just the plain gravitational instability theory, leaving complex non- 
gravitational physics on scales below a Megaparsec to cosmogony, directly in- 
volving star formation. We discuss this possibility below. 
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1.4. What you get is what you see 

An obvious alternative to environmental cosmology was recently discussed by 
Peebles (1999), who pointed out that "as Kuhn has taught us, complex inter- 
pretations of simple phenomena have been known to be precursors of paradigm 
shifts" and perhaps after fifteen years of attempts to salvage the CDM model 
with biasing, it is time to abandon this approach, as well as the biasing hypoth- 
esis itself as "another phlogiston" (all quotations in inverted commas are from 
Peebles 1999). Instead, one can explore a simpler option, that galaxies trace the 
mass distribution, or 

i g = £ and R = b = 1 , (3) 

at least for local (low redshift), optically selected galaxies with a broad magni- 
tude sampling. This approach rests on the idea that no matter how or where 
galaxies form, they must eventually fall into the dominant gravitational wells 
and therefore trace the underlying mass distribution (see Peebles 1980, hereafter 
LSS; Fry 1996). Our purpose here is to test this idea, using measurements of 
relative velocities of pairs of galaxies and the shape of the two-point correlation 
function. 

1.5. Outline of the paper 

In this paper we propose two new tests of the biasing hypothesis, which involve 
two measures of clustering. The first is the two-point correlation function of 
mass density fluctuations, £. The second is a measure of the rate of gravita- 
tional clustering - the relative velocity of particle pairs, v 12 . We describe our 
theoretical model in the next section. Our analytic formulae used to test the 
biasing hypothesis are checked against N-body simulations in §3. In §4 we apply 
our tests to the APM galaxy survey. Finally, in §5 we summarize and discuss 
our results. 



2. Theory 

2.1. The relative velocity 

The relative pairwise velocity v 12 was introduced in the context of the BBGKY 
theory (Davis and Peebles 1977), describing the dynamical evolution of a col- 
lection of particles interacting through gravity. In this discrete picture, v 12 is 
defined as the mean value of the peculiar velocity difference of a particle pair at 
separation f (LSS, Eq. 71.4). In the fluid limit, its analogue is the pair-density 
weighted relative velocity [Fisher et al. 199~4 , Juszkiewicz et al. 1998] , 



w _ ^W 2 )(i + ^)(i + ^)) ( 

where v A and 5 A = p A / (p) — 1 are the peculiar velocity and fractional density 
contrast of matter at a point A = 1,2,.... The brackets (...) denote ensemble 
averages for pairs of points at a fixed separation r = jf^— r 2 |, while £(r) = (5i<5 2 ). 
In gravitational instability theory, the magnitude of v 12 {r) is related to £(r) 
through the pair conservation equation (LSS, Eq. 71.6). 
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Recently Juszkiewicz et al. (1999, hereafter JSD) have proposed an approx- 
imate solution of the pair conservation equation. Using Eulerian perturbation 
theory, they solved the equation for v 12 (r) to second order in £. They also pro- 
posed an interpolation between their large-r perturbative limit, and the well 
known small separation limit - the stable clustering solution, v 12 {r) = —Hr, 
where H is the Hubble constant (LSS, §71). The resulting ansatz is given by 

v 12 {r) = -§ff/rf(r)[l + af(r)] , (5) 

£(r) = (3/r 3 ) f x 2 dx = f(r) [ 1 + £(r) ] . (6) 
Jo 

Here a is a parameter, determined by the logarithmic slope of while 
/ = d InD/dlna, with D(a) being the standard linear growing mode solu- 
tion and a - the cosmological expansion factor (LSS, §11). The solution (||) 
assumes Gaussian initial conditions, and the dynamics of clustering is assumed 
to be dominated by the gravity of inhomogeneities in a pressure-less fluid of 
non-relativistic particles. For all such models, including those with a non-zero 
cosmological constant, / w f2^ 6 (Peebles 1993, §13). If the correlation function 
is given by a pure power law, £ oc r~ 7 , where < 7 < 3, the parameter a is 
given by 

a « 1.2 - 0.65 7 • (7) 

This formalism can be generalized for a non-power law £(r) by replacing 7 in 
Eq. (0) with an effective slope, 

7o = -d]n£/dlnr\ ro , (8) 

evaluated at separation r = r Q , defined by the condition 

C(r ) = 1 • (9) 

JSD tested the equations (||) - (||) against high resolution N-body simulations, 
provided by the Virgo consortium (Jenkins et al. 1998). They found an excellent 
agreement between the streaming velocity, predicted by their ansatz and the 
v 12 (r), measured from the simulations in the entire dynamical range, 0.1 < 
£ < 10 3 . However, the N-body experiments they used were confined to four 
different CDM-like models, considered by Jenkins et al. (1998). As we have 
already pointed out in the Introduction, models of this kind fail to reproduce the 
observed £ 9 (r) unless one resorts to a highly contrived, scale- and time-dependent 
biasing function (Gaztahaga 1995, Jenkins et al. 1998, Peebles 1999). One of 
our objectives here is to test the validity JSD ansatz for v 12 {r) against a new set 
of N-body simulations, which differ significantly from those originally considered 
by JSD. Here we use simulations with a more realistic £(r), inferred by Baugh & 
Gaztahaga (1996) from the measurements of galaxy-galaxy correlations in the 
APM survey (see the description below; from now on, we will refer to these 
numerical experiments and their initial conditions as APM-like). 

2.2. The inflection point 

In the gravitational instability theory, newly forming mass clumps are generally 
expected to collapse before relaxing to virial equilibrium. If this were so, |f 12 ( r ) I 
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would have to be larger than the Hubble velocity Hr to make v 12 (r) + Hr 
negative. As a consequence, the slope of £, 



dln£(r)/cZlnr = -7(f), (10) 

must increase at separations where £(r, t) ~ 1. This effect was recognized long 
ago by Gott & Rees (1975). When the expected "shoulder" was not found in the 
observed galaxy-galaxy correlation function, Davis & Peebles (1976) introduced 
the so-called pre-virialization conjecture as a way of reducing the size of the 
jump in 7(r) (the conjecture involves non-radial motions within the collapsing 
clump; see the discussion in LSS, §71 and Peebles 1 993, pp. 535 - 541; see 
also Villumsen & Davis 1986; Lokas et al. 1996 and Scoccimarro fc Friemanj 



1996j ), Later observational work showed a clear break in the shape of £ for 



several redshift and angular catalogues, which was early evidence for the linear 
to non-linear transition, pointed out by Guzzo and collaborators (see the review 
by Guzzo 1997). 

The arguments, raised by Peebles and Davis (1976) were qualitative rather 
than quantitative. Quarter a century later the precision of N-body simulations as 
well as the quality of the observational data have improved dramatically enough 
to justify a reexamination of the problem. The actual shape of the correlation 
function near £ = 1 can be investigated with high resolution N-body simulations 
like those run by the Virgo Consortium (Jenkins et al. 1998). In all four of the 
Virgo models, the slope of £(r) exhibits a striking feature. Instead of a shoulder, 
or a simple discontinuity in j(r), however, £(r) has an inflection point, 

d 2 £{r)/dr 2 = , (11) 

which occurs at a uniquely defined separation r = r*. At this separation, the 
logarithmic slope of £ reaches a local maximum, dln^/dlnr = —7*. In all 
four of the models JSD investigated, the inflection point indeed appears near 
the transition £ = 1, as expected by the earlier speculations, involving the 
"shoulder" in £. Namely, r* is almost identical with the scale of nonlinearity: 

n « r . (12) 

More precisely, a comparison of Figure 1 in JSD with Figure 8 in Jenkins et al. 
(1998) gives 

\r Q — r*| < 0.1 r a (13) 

for all four considered models. Moreover, for all models, studied by JSD, the —7 
vs. r dependence can be described as an S-shaped curve, with a maximum at 
r = r* ~ r Q , and a minimum at a smaller separation. For r > r*, the nonlinear 
slope (measured from Virgo simulations) follows the linear 7(7*), determined by 
the initial conditions. The separation r* is therefore also the branching point 
between the linear and nonlinear 7(7-) curves, which are identical for r > r* 
(actually, they differ a little; the small differences between the two curves can 
be entirely attributed to sampling errors in the N-body experiments; see Figure 
1 in JSD). This property can be used to identify r* in noisy simulations, and/or 
observations, when the noise in the measured 7(7") curve does not allow a direct 
determination of r* as the position of the maximum in —7(7"). 
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If the relation (12) is indeed a general property of gravitational clustering, it 
can be used as a test of biasing as follows. Suppose the biasing factor is a scale- 
independent constant, significantly greater than unity: 6 > 1. Then £ g 3> £ 
and the relation (|12| ) will break down. For a power-law correlation function, 
£g( r ) = ( r o/ r ) 7 = b^(r), and instead of equation ( O ) we will have 

r* « r Q b- 2h . (14) 

Since the observed slope is 7 ~ 1.8, for 6 = 2, the shoulder in the correla- 
tion function should appear at a separation smaller than a half of the r G ! The 
above argument can be generalized to a broader class of models, allowing scale- 
dependent bias provided b(r) is a smooth monotonic function and b(r ) is sig- 
nificantly greater than unity. Whenever these conditions apply, strong biasing 
always implies r Q S> r*. Hence, the comparison of these two scales determined 
from observations can be used as a diagnostic for biasing. 



3. N-body simulations 
3.1. Initial conditions 

In this section, we compare our ansatz for v 12 (r) with the results from P 3 M 
simulations. We use APM-like models for the initial shape of P{k) in Baugh 
& Gaztahaga (1996). The box size is 600 h' 1 Mpc (or 300 h' 1 Mpc) with 
200 3 (or 100 3 ) dark matter particles. The APM-like models have Gaussian 
initial conditions with a power spectrum inferred from correlations in the APM 
survey, following the procedure, introduced by Baugh & Gaztahaga 1996. The 
power spectrum of the density fluctuations is related to £(r) through the usual 
Fourier transform, 

P(jfe) = J £(r)exp(ik-r)d 3 r . (15) 

Following the prescription of Baugh & Gaztahaga, we assume an initial power 
spectrum of the form 

Ak 

P(k) = re • (16) 

\ + (k/0mh Mpc" 1 ^ 



The linear spectrum, given above is designed to evolve into a nonlinear one, 
matching the APM measurements of P{k) under the assumption of no bias. The 
normalization constant A is directly related to the degree of inhomogeneity of 
the mass distribution at the end of the simulation, parametrized by cr 8 - the rms 
mass density contrast, measured in spheres of a radius of 8 /i -1 Mpc. Following 
Baugh and Gaztahaga, we choose the constant A in order to ensure that at the 
end of the simulation, a 8 = 0.85. The quality of the agreement of the evolved 
£(r) with the APM data at small separations, where £ S> 1, depends on the 
assumed value of O m . However, as we show below, for r > 2/i _1 Mpc and £ < 3, 
this effect is negligible (compare Figures 1 and 3). This range of separations and 
clustering amplitudes brackets from both sides the region, on which we focus on 
here: the transition between £ < 1 and £ > 1. We use simulations with A = 
and two different values of the density parameter: Q m = 1 and 0.3. 
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3.2. Estimators 



Two different estimators for £(r) and v 12 are used. To construct the first estima- 
tor, we cut the simulation box into cubical pixels of size Aj, placed on a regular 
lattice. The density fluctuation amplitude at the ith pixel is 

X 1, (17) 



(TV) 

where Ni is the particle count in that pixel. The estimator for the two-point 
function is then: 

= jEWiW > ( 18 ) 

r i,j 

where 

N r = £Wi,-(r) (19) 

i,j 

is the number of pairs of pixels at separation r in the sampled region, and the 
window function Wij (r) = 1 if pixels i and j are separated by |r*j — fj \ = r ± Ar , 
and otherwise. For the pairwise velocity we define as V{ the average velocity 
in the ith pixel at position fl. We can then use an equivalent expression for the 
first estimator of v 12 : 

t)U(r) = ^(vi-v^-rijUijir), (20) 



where 



(21) 



is a unit vector separating pixels i and j and the sum is over all pairs of pixels, 
while 



_^ (1 + 5,) (l + ^W^r) 

N r [i+ary 



The second estimator uses all of the particle pairs rather than the counts in pixel 
pairs: 

i (2) (r) = t4fE WUt) -1, (23) 



(n)V r 



where now the window function W(j(r) = 1 if particles i and j are separated 
by \Ti — fj\ = r ± Ar, (n) is the mean particle density and V r is the volume of 
the spherical shell of radius r and thickness Ar. Typically, (n) V r is estimated 
from a random catalogue of particles, drawn from the same sample, because 
the effective volume, containing the pairs separated by a distance r ± Ar might 
depend on the geometry. However, our simulations use a large and periodic box, 
with no boundaries and high densities (not affected by shot-noise at the scales 
of interest). So the above expression gives a good and quick estimator. The 
corresponding estimator for the pairwise velocity is 

9 



where V{ and Vj are now individual particle velocities. Note that this estimator 
does not depend on the effective volume of the shell, V r . 

Both estimators agree reasonably well in our simulations. The first set is 
more useful (faster to run) for large separations, as we can reduce the resolution 
of the lattice and have a relatively small number of pixels. The second set is 
more adequate (faster to run) for the small separations. 

3.3. The correlation function 

The evolved, nonlinear correlation functions, measured from simulations are 
shown in Figure |l| (top panel). The full squares correspond to the fL m = 1 
model, while the open squares represent £l m = 0.3. For comparison, we show 
the linear correlation function (dashed line), scaled from some "initial" redshift 
z Q to the present (z = 0), following the linear theory expression for the Einstein- 
de Sitter model, £(r, 0) = £(r, z ){l+z ) 2 . Nonlinear effects are more pronounced 
in the low density model. This happens because of the well known suppression 
of linear growth, which occurs at late times (z < l/O m ) in low density models 
and leads to the enhanced clustering on small scales relative to large scales. 

Note however, that although the correlation functions differ significantly in 
amplitude at separations r < 2/i _1 Mpc, their slopes 7(7") are almost indistin- 
guishable (Figure [l], bottom panel). 

The particle resolution (the Nyquist wavelength cx iV~ 1//3 ) of the simulations 
used here is significantly lower than the resolution of Virgo simulations, and the 
noise in the measured £(r) is further amplified by differentiating over r. As a 
result, determining the position of the inflection point r* directly from the 7(7") 
curve alone is difficult. To overcome this problem, we identify r* by comparing 
the linear and nonlinear 7(7") curves. Taking r* to be the separation at which 
the nonlinear slope drops below the linear slope in Figure [| as described earlier, 
we get 

r* ~ 5 /T^Mpc sa r a , (25) 

in excellent agreement with equation (|l~3|). Hence, the equality between r* and 
r can probably be considered as a generic outcome of gravitational dynami- 
cal evolution in a model where galaxies trace the mass and the initial slope, 
cf ln£/dlnr, is a smooth decreasing function of the separation r. Such a picture 
is also known as hierarchical clustering; an obvious additional condition to make 
sure that small scale clumps collapse before the large scale ones, is 7 > 0, see 
e.g. LSS. 

3.4. The relative velocity 

In this section we describe N-body tests of the JSD model for the relative 
motions in pairs of galaxies. We consider two models with APM-like initial 
spectra: an Einstein-de Sitter model and an open model with £l m = 0.3. For 
both models, the theoretical predictions for the mean pairwise velocity, based 
on equation (|5|), are plotted in Figure [2] as continuous lines. These predictions 
can be compared with N-body measurements, shown as full squares for the 
Q m = 1 model and as open squares for Q m = 0.3. The agreement between 
the theoretical and experimental v 12 (r) curves shows that our ansatz provides a 
good approximation of the N-body results in the entire dynamical range probed 
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Figure 1. The top panel shows the linear correlation function (dashed 
line) and the measured non- linear £(r), obtained from the APM-like 
simulations with density parameters Q m = 0.3 (open squares) and 
VL m = 1.0 (full squares). The bottom panel shows the correspond- 
ing logarithmic slope, —7(7") = dln^/dlnr for each of the three curves 
from the top panel. The vertical dotted line shows the separation r Q , 
defined by the condition £(r ) = 1 (top) and the separation r*, at which 
the non-linear j(r) curve crosses the linear one (bottom). 
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Figure 2. The mean pairwise velocity v 12 (r), measured from two sets 
of APM-like simulations with O m = 0.3 (open squares) and Q m = 
1.0 (full squares), are compared with with an approximate analytical 
solution of the pair conservation equation (eq. ||; continuous lines). 

for both models. The mean and errors in the O m = 1 simulations come from 
the mean and dispersion, obtained from five independent realizations of the 
APM-like model. For the open model (ft m = 0.3, open squares) we used only 
one realization, but the expected sampling variance is expected to be the same. 
Indeed, the initial P(k) is identical in both cases. Moreover, the long-wave tails 
of the final power spectra (which determine the size of the sampling error bars) 
are also identical because they are not affected by the nonlinear evolution. 

4. Comparison with observations 
4.1. The correlation function 

The measurements of Cg( r )i an d j{r) = — <iln£ 9 /<ilnr, obtained from the angu- 
lar correlations of galaxy pairs in the APM catalogue (Baugh 1996), are plotted 
in Figure ^. The top panel shows the two-point function (points with error bars), 
and the linear theory curve, described in §3.3 (dashed line). The intersection of 
the two perpendicular dotted lines marks the point (£ 9 , r) = (1, r Q ). The bottom 
panel of Figure ||| shows the APM 7(7-) as a function of the pair separation r. 
Note the remarkable similarity between the empirical data and the characteristic 
peak in the j(r) found in the simulations (Figure |l|). The intersection of the 
two mutually perpendicular, dotted lines in the bottom panel of Figure ^ marks 
the crossing between the linear model for 7(7") (dashed line) and the nonlinear 
7(7") curve, determined from the APM catalogue. The crossing occurs at the 
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Figure 3. The spatial correlation function of APM galaxies (top 
panel, symbols with errorbars), compared to the linear theory APM- 
like model, described earlier (top panel, dashed line). The bottom panel 
shows the corresponding logarithmic slope, 7(7*) . The intersection of 
the two perpendicular dotted lines marks the points where £ ~ 1 (top) 
and where the non-linear slope crosses the linear one (bottom). 
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separation r ~ 5h~ Mpc, and to first approximation this scale could be identi- 
fied with r*. However, a closer inspection of our Figure || suggests that, given 
the error bars, the actual position the peak could be shifted to the right, to a 
somewhat larger separation. Taking into account the error bars in Figure [3| as 
well as the uncertainties in the assumed linear theory slope, we obtain 

r* ~ (6±1) h^Mpc , (26) 

and 

r Q ~ (5±1) /» -1 Mpc. (27) 

The slope at r = r* is 7* ~ —1.4. If we assume the linear bias model, the 
relation ~ r 6~ 2 / 7 * gives 

b = 1.11 ± 0.22 (28) 
at one-sigma statistical significance level. 
4.2. The relative velocity 

We will now apply the second of the two proposed tests of biasing: the relative 
velocity test. We will compare the mean pairwise velocity, predicted by assuming 
that the APM galaxies trace the mass with the pairwise velocity, measured 
directly from a peculiar velocity - distance survey. 

Figure H shows v 12 (r) curves, predicted by equation (|J) (continuous lines) for 
three different values of O m , from bottom to top U m = 1.0,0.3,0.1, respectively. 
To calculate v 12 (r), we have used £(r), estimated from the APM survey under 
the assumptions £ s = £ and R = 1. 

Before making the comparison, we must overcome the following problem. 
The survey has a significant depth, with the mean redshift of z ~ 0.15 while 
the measured v 12 (r) corresponds to the present time [z = 0). To evolve £(r, 2) 
from z ~ 0.15 to z = 0, we need to make some additional model assumptions. 
Gaztahaga (1995) has shown that for this redshift range, the uncertainties in 
the details of dynamical evolution of £ are small. In particular, choosing an 
incorrect value for Q m can affect £ at most at the several per cent level (for 
O m ranging from 1 to 0). Adding this to other possible sources of errors, such 
as uncertainties regarding the redshift evolution of the galaxy number density 
and sampling and selection fluctuations, Gaztahaga (1995) estimates that the 

resulting relative uncertainty in the amplitude of £ g is ~ 20%. According to 
his analysis, the present {z = 0) amplitude of the rms fluctuation of the APM 

galaxy counts, measured in spheres of radius of 8 /i -1 Mpc, is 1.1 ~ a s APM ~ 0.9. 

To be conservative, for each value of f2 m , we plot the predicted v 12 (r) curves 
for two values of <7 g , differing by 20%. The resulting prediction for each value 
of 0, m is therefore an area rather than a single v 12 (r) curve (see Figure 4). The 
lower boundary of each shaded area assumes a s = 1.1 while the upper boundary 
was calculated by assuming cr 8 = 0.9. A direct measurement at r = 10 h^ 1 
from the Mark III galaxy peculiar velocity survey (Willick et al. 1997) gives 
KJuszkiewicz et al. 2000a) 

v 12 = -280 ± 60 km/s . (29) 
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Figure 4. Predictions for the mean pairwise velocity v±2, based on 
the assumption that the APM galaxies trace the mass. The shaded 
regions correspond to 20% uncertainties in the strength of clustering. 
We consider three values of O m = 1.0,0.3,0.1, top to bottom (as la- 
beled). The point with error bars corresponds to a direct measurement 



from the Mark III survey (Juszkiewicz et al. 2000a) 



It is reassuring that this value, plotted in Figure 4, overlaps with the shaded 
area, corresponding to O m = 0.3, because it agrees with ranges for fi m and a s 
obtained by Juszkiewicz et al. from the analysis of the Mark III survey alone. 
Their one-sigma constraints are VL m = 0.35i *25 an d °8 > 0.7, and the analysis 
assumes that the correlation function for the mass is well approximated by a 
pure power law with 7 = 1.75. 

From the agreement between the predicted and observed value of the mean 
pairwise velocity, we conclude that the Mark III and APM data, taken together, 
are consistent with the hypothesis that the APM galaxies trace the mass, b ~ 
R ~ 1, while the density parameter is low, Q m ~ 0.3. 

4.3. How can biasing affect v 12 (r) ? 

Apart from leading to predictions which are confirmed observationally, the "what 
you get is what you see" hypothesis has another important advantage: simplicity. 
Once £(r) is estimated from observations, a family of v 12 (r, Q m ) curves can be 
calculated directly from equation (|5|) for any given range of values of the density 
parameter. This simplicity will immediately go away if we allow scale-dependent, 
stochastic and nonlinear biasing. A frank answer to the question posed in the 
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heading of this subsection would have to be "Only God (of biasing) knows". 
Predicting v 12 {r) without resorting to massive numerical simulations would be 
simply impossible. We can get the idea of what is in store by considering only the 
leading order term in the perturbative expansion for v 12 (r) at large separations 
(Juszkiewicz et al. 2000b), 



v 12 (r) = -lHrf(Q m )i gp (r) , (30) 

where £, gp {r) = (3/r 3 ) J Q r £, gp (x)x 2 dx is the galaxy-mass cross-correlation func- 
tion, 

C ap (r) = (5(0)5 g (f)) , (31) 

averaged over a sphere of radius r. The function £ gp describes the cross-correlations 
between the mass density and the density of galaxies in the velocity field survey, 
which in the case considered here would be the Mark III catalogue. 

To make progress in our analysis, we will now generalize the definition of 
the stochasticity parameter introduced as a normalized cross-correlation of two 
random fields, 5 and 5 g , measured at the same position in space. Instead, we 
will consider a cross-correlation of the same two fields measured at two different 
positions in space, separated by distance r. Our old equation (||) is replaced by 

#( r ) = JkM . (32) 



Let us make another simplifying assumption, that b 2 (r) = £ 5 (r)/£(r), as well as 
R, are separation-independent. Equation (|30|) becomes 



vy 2 (r) = -lf(n m )HrRb^r) . (33) 

The expected relative pairwise velocity can now be related to the APM data by 
substituting 

= FTTT f U^dx , (34) 
b A r 6 Jo 

where b A 2 = £ s (?*)/£( r ) and £ g is the APM galaxy correlation function. In case 
of trouble in predicting the correct v 12 {r) curve, we now have three essentially 
free parameters which can be readjusted. This is only the tip of the iceberg, as 
we have ignored nonlinear dynamics as well as the scale-dependence of b and R. 

For the linear bias model, the predictions are in clear conflict with the data 
unless b is close to unity. After setting R = 1, we get v 12 (r) oc b. Then, if the 
biasing factor for spiral galaxies is, as usually assumed b w 1, our predictions for 

v 12 (r) in the linear regime (r ~ 10/i _1 Mpc) should be similar to the unbiased 
predictions already plotted in Figure (Q) . If the biasing factor for the ellipticals 
is significantly different, say, b ~ 2, the elliptical subsample of the Mark III 
data should give estimates of v 12 which differ from the estimates from the spiral 
sample by the same factor of two. Meanwhile the estimates from the appropriate 



subsamples in the real data are indistinguishable (Juszkiewicz et al. 2000a) 
Hence, just as in case of the shape of the APM correlation function, considered 
above, the deterministic linear biasing model is inconsistent with observations. 
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We can summarize the last two subsections as follows. The prediction for 
v 12 (r), based on the assumption that the APM galaxies trace the mass passes 
our test as it agrees with the velocity, estimated from the Mark III data. The 
simplest prescription of biasing fails the test. More complicated prescriptions 
can probably be made to pass, which is not surprising, given the number of free 
parameters. 

5. Summary and discussion 

A quarter of a century ago, Gott and Rees predicted that gravity should leave 
its mark on the shape of the galaxy autocorrelation function: a "shoulder" , or 
steepening of the slope of the correlation function should appear near the sep- 
aration where £ passes through unity. At the time, biasing was unheard of, 
and Gott and Rees (1975) assumed £ = £ 9 . Recently, in another context, JSD 
have studied the the £ = 1 boundary in the evolution of the mass correlation 
function, using results from Virgo simulations. They found that the "shoulder" 
is actually an inflection point, occurring at a well defined separation r*. In all 
four CDM-like models they studied, the nonlinear transition looked strikingly 
similar: the inflection occurred at almost the same separation as that of the 
nonlinear transition: r* ~ r c , where r Q corresponds to £ = 1. Here we have 
tested the degree of universality of their results by widening the range of mod- 
els considered. Our additional objective was to study the range of validity of 
an approximate solution of the pair conservation equation, proposed by JSD 
to study the nonlinear evolution of the relative velocity of particle pairs at a 
fixed separation, v 12 (r). We used N-body simulations, with APM-like initial 
conditions, with two different values of the density parameter: Q m = 1 and 0.3. 
The APM-like initial power spectra differ significantly from all of the CDM- 
like spectra, considered earlier by JSD. Moreover the spectra of the latter kind 
appear as more realistic to us because they can reproduce observations without 
resorting to scale-dependent biasing. Our APM-like simulations are in excel- 
lent agreement with earlier results, confirming the validity of the JSD ansatz for 
v 12(f) and the conjecture that the appearance of the shoulder in the correlation 
function near the £ = 1 transition is a feature of gravitational dynamics rather 
than a peculiarity of a particular set of initial conditions. 

Using these results, we proposed two tests of the hypothesis that galaxies 
trace mass. The first of the tests is based on an obvious idea, that if £(r) = 
£g(r), the galaxy correlation function near £ ff = 1 should exhibit properties 
similar to those of the matter correlation function. We examined the behavior 
of the correlation function, derived from the APM catalogue and found exactly 
the same features we knew earlier from N-body simulations, in particular the 
agreement between the two characteristic scales, w r a . It is difficult to imagine 
how such an agreement could happen by a mere coincidence, which would have to 
be the case if £ ff is unrelated to £. The agreement between the two characteristic 
scales can be used to constrain the linear biasing factor for the APM catalogue 
to be within 20% of unity. This constraint agrees with an earlier limit, obtained 
from measurements of the three-point correlation function from the APM survey 
(Gaztanaga 1994, Frieman & Gaztahaga 1999). 
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The second test confronts the v 12 , predicted by assuming that the APM 
galaxies are unbiased tracers of mass with direct measurements of v 12 . The 
results are again consistent with b ~ 1 and a low density parameter, £l m ~ 0.3, 



in agreement with the limits obtained from the velocity data alone (Juszkiewicz 



et al. 2000a 



We are impressed how well the observations are reproduced by the simple 
calculations based on the assumption that galaxies follow the mass distribution, 
at least on large (weakly non-linear) scales. We are unable to constrain biasing 
models with a large number of free parameters, but their predictive power is 
questionable and one may ask: are such models falsifiable and therefore worth 
constraining? 

Our results are by no means final, they are also less rigorous than one could 
wish because we are limited by the accuracy of the present observational data. 
New generation of catalogues promise a dramatic improvement on this front in 
the near future (for an excellent collection of reports on the state of the art in 
this field, see Colombi et al. 1998). 
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